Abstract. We present the results of a detailed non- 
LTE analysis of the ultraviolet and optical spectrum 
of the 06.5 Iaf+ star HD 153919 - the mass donor 
in the high-mass X-ray binary 4U1700-37. We find 
that the star has a luminosity log(L*/Lo)=5.82±0.07, 
Teff=35,000±1000 K, radius R,=21.9lJ;^ Ro, mass-loss 
rate M=9.5xl0~^ Mq yr~^, and a significant overabun- 
dance of nitrogen (and possibly carbon) relative to solar 
values. Given the eclipsing nature of the system these re- 
sults allow us to determine the most likely masses of both 
' components of the binary via Monte Carlo simulations. 
O ; These suggest a mass for HD 153919 of = 58 ± 11 Mq 
- implying the initial mass of the companion was rather 
high (^60 M0). The most likely mass for the compact 
' companion is found to be M^; =2.44±0.27 Mq, with only 
' 3.5 per cent of the trials resulting in a mass less than 
\^ ' 2.0 Mq and none less than 1.65 Mq. Such a value is sig- 
nificantly in excess of the upper observational limit to the 
masses of neutron stars of 1.45 Mq found by Thorsett & 
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Chakrabarthy (1999), although a mass of 1.86 Mq has 
' recently been reported for the Vela X-1 pulsar (Barziv 
' et al. [200 1| ). Our observational data is inconsistent with 
CO ' the canonical neutron star mass and the lowest black hole 
mass observed (^4.4 Mq; Nova Vel). Significantly chang- 
' ing observational parameters can force the compact object 
I mass into either of these regimes but, given the strong pro- 
portionality between and M^: , the 0-star mass changes 
_ , by factors of greater than 2, well beyond the limits deter- 
I ■ mined from its evolutionary state and surface gravity. The 
O I low mass of the compact object implies that it is difficult 
• to form high mass black holes through both the Case A 
^ ■ & B mass transfer channels and, if the compact object 
• • ■ is a neutron star, would significantly constrain the high 



>< 



density nuclear equation of state. 
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1. Introduction 



First detected by the Uhuru satellite (Jones et al. 1973) 
the eclipsing X-ray source 4U 1700-37 was quickly associ- 
ated with the luminous 06.5 Iaf+ star HD 153919, con- 
firming 4U 1700-37 as a high-mass X-ray binary (hence- 
forth HMXB). HMXBs are systems composed of an OB 
star and a compact companion (neutron star or black 
hole), with the X-ray emission resulting from the accretion 
of material by the compact companion. In the subclass of 
supergiant HMXB systems material is accreted either via 
Roche lobe overflow or directly from the powerful stellar 
wind of the OB primary. Given that HD 153919 slightly 
underfills its Roche Lobe (e.g. Conti 1978 ) mass transfer 
proceeds via the latter mechanism. 

Alth ough an orbital period of ~3.412 days (Jones et 
al. 1973 ) was quickly identi fied fo r 4U 1700-37, extensive 
searches (e.g. Rubin et al. 1996 and references therein) 
have failed to identify any other X-ray periodicities within 
the system that might correspond to the pulse perio d for 
a possible neutron star (although Konig & Maisack 1997 
claim the presence of a 13.81 day period in CGRO BATSE 
& RXTE ASM datasets). Given the absence of any X-ray 
pulsations and the unusually hard nature of the spectrum 
various authors (e.g. Brown, Weingartner & Wijers 1996) 
have suggested that the compact companion could be a 
low mass black hole rather than a neutron star. However, 



Reynolds et al. ( |1999| ) point out that the 2-200 keV spec- 
trum of 4U 1700-37 differs from those commonly observed 
for black hole candidates such as Cygnus X-1. Given that 
the X-ray spectrum of 4U 1700-37 is qualitatively similar 
to those of accreting neutron stars they suggest that the 
compact object is also a neutron star, and explain the lack 



of pulsations as due to either a weak magnetic field or an 
alignment of the magnetic field with the spin axis. 

With a spectral type of 06.5 Iaf+, HD 153919 is the 
hottest and potentially most massive mass donor of any 
of the the HMXB systems. As such, determination of its 
fundamental parameters (radius, temperature, mass and 
chemical composition) is of importance given that these 
will potentially provide valuable insights into the evolution 
of very massive stars and their ultimate fate. In particular 
by determining the masses of the components of HMXB 
systems, limits to the progenitor masses for neutron stars 
and black holes in such systems can be found|^ while the 
chemical composition of the mass donor can shed light on 
the pre-supernova (SN) mass transfer mechanisms. 

The paper is ordered as follows. Section 2 describes the 
determination of the physical properties of HD 153919; the 
complete dataset and the non-Local Thermal Equilibrium 
(non-LTE) code used to analyse it. Section 3 describes 
the Monte-Carlo technique used to determine the masses 
of both components of the binary system. In Sections 4 & 
5 we discuss the implications of our results for the evolu- 
tion of hot massive stars, limits for the progenitor masses 
of compact objects and the equation of state for nuclear 
matter. Finally in Section 6 we summarize the main re- 
sults of the paper. 



2. Determination of the stellar parameters for 
HD 153919 

Determination of the fundamental stellar parameters of 
HD 153919 is complicated by its high temperature and 
mass loss rate, which necessitates a sophisticated non- 
LTE treatment. Early attempts to determine the physi- 
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* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (64.H-0224) 



we note that Wellstein & Langer ( 1999 ) demonstrate that 
such limits derived from binary systems cannot be directly ap- 
plied to single stars 
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cal properties of HD 153919 suggested that the star might 
be undermassive by a factor of 2 (e.g. Conti 1978, Hutch- 
ings 1974) which would imply that the terminal velocity 
of the stellar wind (Voo ~1700 km s~^; van Loon, Kaper 
& Hammerschlag-Hensberge 2001 (henceforth vLOl)) is a 
factor of 10 in excess of the escape velocity (Vcsc)- Given 
that Howarth & Prinja ( |1997D show that Voo/Vesc < 4 
for O stars (and typically ~2.5) this discrepancy clearly 
needs to be resolved. More recent analys es sti ll do not re- 
solve the issue, with Heap & Corcoran (|l992| ) suggesting 
a mass for HD 153919 of M*=52±2 Mq (thus broadly in 
line wi th th e expected mass for such a star) while Rubin 
et al. ( |1996D propose M,=30t" Mq, suggesting that the 
star is probably undermassive for its spectral type. 

As will be shown in Sect. 3, determination of the 
masses of both components in the system is hampered 
by the considerable uncertainties in the stellar radius of 
HD 153919. In order to address this problem, and in the 
light of dramatic advances in the sophistication of non- 
LTE model atmospheres we have decided to reanalyse 
both new and published ultraviolet to near-infrared spec- 
troscopic and optical to mid-infrared photometric obser- 
vations of HD 153919 in order to refine previous estimates 
of the stellar parameters. 



2.1. The complete dataset 

Archival and new ultraviolet to mid-infrared spectroscopic 
and photometric data were used to derive a set of stellar 
parameters for HD 153919. High spectral resolution UV 
spectroscopy was obtained with the International Ultra- 
violet Explorer (lUE); the data used and reduc tion p roce- 
dures employed are described in Kaper et al. (1993); the 
details are not repeated here. Four high S/N and high 
spectral resolution (R=48,000) optical spectra (~3700- 
86OOA) were obtained in 1999 April with the Fiber-fed 
Extended Range Optical Spectrograph (FEROS) mounted 
on the ESO 1.52m telescope at La Silla. All 4 spectra 
were wavelength calibrated and optimally extracted to de- 
termine if significant changes in the spectrum occured at 
different orbital periods. Besides line-profile variability in 
the strongest "wind" lines and the shift in radial veloc- 
ity due to orbital motion, no evidence is found for in- 
trinsic variability of the photospheric spectrum. The fi- 
nal spectrum used for determining the stellar parame- 
ters of HD 153919 was that taken during X-ray eclipse 
to further minimize the effects of any perturbation of 
the wind by the presence of the compact companion 
(see Sect. 2.2). Near- infrared spectra between 1-2.2/^m 
and optical to near-infrared photometry were taken from 
Bohannan & Crowther ( 1999| ) and mid-IR photometry 
( 6.8/im =669 mJy, 11.5/im=244 mJy) from Kaper et al. 
(1997); see respective papers for the particular reduction 
strategies employed in each case. 



Fig. 2. Plot of the observed (solid line and data points) 
and theoretical (dotted line) UV - mid-IR spectral energy 
distribution for HD 153919 




% 



HD 153919 
E(B-V) = 0.53 



3.00 3.25 3.50 3.75 4.00 4.25 4.50 4.75 5,00 

log Wavelength (A) 



2.2. Spectral analysis 

To determine the stellar properties of HD 153919 we have 
utilised the non-LTE code of Hillier & Miller ( |1998D which 
solves the radiative transfer equation subject to the con- 
straints of statistical and radiative equilibrium, in a spher- 
ical, extended atmosphere. Line blanketing is incorporated 
directly through the use of a super-level approach. We use 
a similar atomic model to that employed by Crowther et 
al. ( ^002 ) in their study of early O supergiants, including 
Hi, Hei-ii, Ciii-iv, Niii-v, Oiii-vi, Siiv, Piv-v, Siv-vi 
and Feiv-vii. For extreme O supergiants, line blanketing 
and the strong stellar wind conspire to produce significant 
differences in stellar parameters relative to the sta ndard 
plane-parallel hydrostatic results (see Crowther et al. 2002 
for further details). 
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Our procedure is as follows. We adjust the stellar tem- 
peratureQ and mass-loss rate of an individual model until 
the 'photospheric' Hen A4542 and Hei A4471 Hues are 
reproduced. Simultaneously, we vary the total mass-loss 
rate until Ha is also matched. The exponent of the /3-law 
is adjusted until the shape of Ha is well reproduced - 
for HD 153919 we obtain /3 ~1.3. The input atmospheric 
structure, connecting the spherically extended hydrostatic 
layers to the /3-law wind is achieved via a parameterized 



scale height, h (see Hillier et al. 2002 for details), for which 
/i=0.001 yields a reasonable match to Hei and Balmer line 
wings, consistent with log 5=3.45-3.55. We adopt a termi- 
nal wind velocity of Woo=175G kms~^ (vLOl; Howarth et 



al. |1997D . 

The formal solution of the radiative transfer equation 
yielding the final emergent spectrum is computed sepa- 
rately, and includes standard Stark broadening tables for 
Hi, Hei-ii. Except where noted, these calculations assume 



^ Defined, as is usual for an extended atmosphere, as the ef- 
fective temperature corresponding to the radius at a Rosseland 
optical depth of 20 



a microturbulent velocity Wturb==10 kms^^. Hillier et al. 
( 2002 ) also discuss the effect of varying Wturb in O star 
models. Additionally, we find good agreement with ob- 
servations using w sin 1=150 kms~^ (Howarth et al. 1997 
derived 120 kms^^). 

It is extremely difficult to determine accurate He/H 
abundances in O supergiants as discussed by Hillier et 
al. ( 2002 ). Consequently, we adopt He/H=0.2 by number, 
whilst C and N abundances are varied until diagnostic 
optical line profiles are reproduced. In Fig. 1 we present 
selected optical line profile fits to FEROS observations 
of HD 153919. Overall, agreement for rcff=35kK is very 
good, with the exception of Hen A4686. Hei A4471 pro- 
vides our main temperature constraint since other blue 
optical Hei lines are weak or absent. Alternatively, we 
considered using Hei A5876 (or A10830) together with 
the Hen A4686 line. However, this method (followed by 
Crowther & Bohannan 1997) yields significantly ('^4kK) 



lower stellar temperatures, and suffers from inconsisten- 
cies involving the ionization balance of UV/optical metal 
lines. Therefore, we have greater confidence in our adopted 
diagnostics, which do not suffer from such problems. 
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Fig. 3. Plots of selected regions of the UV spectrum of HD 153919 (solid line) and best fit model (dotted line); 
parameters listed in Table 1 . Note that our model does not take into account the observed Raman-scattered emission 
lines (which are not of a photospheric origin) in the range 1400-1700 A. 
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From spectral energy distribution fits to lUE spec- 
trophotometry and Johnson photometry (Fig. 2), we de- 
rive Eb-v = 0.53±0.02 mag. Alternatively, using the in- 
trinsic colour of {B — V)() = —0.30 from the early O 
supergiant calibration of Schmidt-Kaler (1982), we de- 
rive Sb-v=0-55 from Johnson photometry of HD 153919 
(V=6.54 and B - V=0.25, Bolton & Herbst |l976|). Con- 
sequently, we adopt i?B-v=0.54±0.02 for the remainder 
of this work. We adopt a d istanc e modulus of 11.4 mag 
to HD 153919 (Ankay et al. ^00l|) , implying My = -6.53 
mag. Our derived temperature provides a bolometric cor- 
rection of —3.3 mag, therefore we obtain log(L/L0)=5.82 
for HD 153919, and thus 7?=21.9 Rq. The derived mass- 
loss rate is 9.5x 10~^ M0 yr~^, assuming that the Ha line- 
forming region is not clumped. Moderate clumping would 
reduce this value by a factor of ~2. Derived parameters 
are listed in Table 1, and are in reasonable agreement with 
those derived by vLOl through the analysis of ultraviolet 
resonance lines based on the Sobolev with Exact Integra- 
tion (SEI) method. 

Our primary nitrogen abundance diagnostics are N ill 
A4634-41 and A4097, which together imply eN — 9eN,0. 
With this value, the very weak N ill A5320-4 feature is well 
matched, but other lines in the vicinity of He 11 A4542 are 
somewhat too strong (likely due to an incomplete treat- 
ment of Niii quartet states in our models). Carbon is 
somewhat more difficult to constrain, with Cm A4647- 
51 weh matched for ec = l-Oec,©- Civ A5801-12 is weh 
reproduced with this value, whilst C ill A5696 is too weak, 
implying a yet hi gher a bundance. As discussed elsewhere 
Hillier et al. 



strongly predicted nor observed in the A1200-1400 region 
(see Crowther et al. 2002 for further details). 

While it is possible to explain the nitrogen enrichment 
in terms of rotational mixing it is impossible to produce 
carbon enrichment via this mechanism. Any carbon pro- 
duced in the helium burning layers of the star has to pass 
through the hydrogen burning layers before reaching the 
surface where it will be converted to nitrogen. Therefore, 
any excess carbon in HD 153919 is therefore likely to result 
from mass transfer from the more evolved binary compo- 
nent prior to SN - this will be returned to in Sect. 4. 

Given the presence of a compact companion for 
HD 153919, it is reasonable to ask whether the assumption 
of spherical geometry is justified - does the X-ray flux lead 
to significant departures from spherical symmetry for the 
ionisation of the wind (which in turn could lead to mod- 
ifications in the line driving force)? Hatchett & McCray 
(1977) suggest that the X-ray emission will lead to a re- 
duction of moderately ionised atoms in the wind (such 
as Si IV and Civ). Given that the ionised zone will move 
with the compact object we might expect to see orbital 
modulation in some of the wind UV resonance lines as the 
ionised zone passes in front and behind the stellar disc (or 
from the presence of a photo-ionization wake in the sys- 
tem, cf. Kaper et al. 1994). However, there is no convinc- 
ing evidence for orbital modulation in the UV resonance 
line due to the Hatchett-McCray effect (e.g. Kaper et al. 



1990, 1993); the small changes in line profiles with orbital 



phase are instead most likely due to Raman scattering of 
EUV photons generated by the X-ray source (Kaper et al. 



1998D, oxygen is exceedingly difficult to Fq^ [[qq%_ Additionally, the modeling was performed on 



(e.g. 

constrain in mid-0 supcrgiants due to lack of suitable op- 
tical diagnostics. The high nitrogen overcnrichment is not 
easily explained via single star evolution, unless carbon 
(and to a lesser degree oxygen) is very depleted via the 
CN (or ON) cycle. Crowther et al. (2002) discuss similar 
problems for Magellanic Cloud O supcrgiants. 

Turning to UV comparisons, we show rectified high 
resolution lUE spectroscopy of HD 153919 (phase 0.15) 
in Fig. 3 together with synthetic spectra. Overall, agree- 
ment for He II A1640, C IV A1550 and N iv A1718 is reason- 
able, with predicted Si IV A1393-1402 emission too weak 
adopting egi — 1-Oesi,0- Since X-rays are not explicitly 
considered in this study, the shocked UV Nv A1238-42 
resonance doublet is predicted to be too weak. The sole 
prominent oxygen feature present in the UV (or optical) 
region is O iv A 1338-43, which is reasonably well matched 
with Co = 0.5eo,o, although we do not claim that this 
represents an accurate constraint. 

Additional, powerful evidence in favour of our derived 
temperature is the good match between the synthetic 
Feiv-v spectrum and observations, again with epe = 
LOepo,©)- Fig. 2 shows good agreement with the dominant 
Fev 'forest' observed between A1300-1600 in HD 153919, 
plus the weaker Fe iv in the A1500-1800 region. Fe vi is not 



the spectrum obtained during the X-ray eclipse to further 
minimise any possible effects of irradiation on the stellar 
wind (cf. Sect. 2.1). 

Using a modified 2-dimensional Sobolev Exact Integra- 
tion (SEI) code vLOl analyse the UV line variability and 
confim that any Stromgren sphere caused by the presence 
of the X-ray source is rather small, and will have a neg- 
ligible effect on the ionization structure and line driving 
of the wind (since the wind is dense and the ionizing fiux 
low). Equally, the Stromgren zone does not extend to the 
surface of the star and so should not lead to a significant 
degree of X-ray heating of the stellar surface. 

Phase resolved continuum observations (e.g. vLOl, 
Hammerschlag-Hensberge & Zuiderwijk 1977, van 
Paradijs et al. 1978) constrain orbital variability to 



<4 per cent in the UV and 4-8 per cent in the optical, 
indicating that HD 153919 shows little departure from 
sphericity (possibly as a result of a large mass ratio). 
Note that continuum emission from the stellar wind is 
essentially negligible at wavelengths shorter than a few 
microns. Therefore, the lack of significant variability 
cannot be attributed to emission from the outer regions 
of the stellar wind 'shielding' a heavily perturbed stellar 
surface and/or inner wind from view. 
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Parameter 



Value 



E(B-V) 

log(L/LQ) 

R. 

M 



0.54±0.02 
35,000±1000K 
5.82±0.07 



21.r 



^1-3 R, 
9.5x10"'^ Mq yr"^ 
1750 km s"^ 
3.45-3.55 



Table 1. Stellar parameters for HD 153919 derived from 
the NLTE modeling described in Sect. 2.2. 



The photometric variability further constrains any 
change in stellar temperature due to X-ray heating to less 
than the uncertainty in the stellar temperature derived 
from our NLTE modeling. Therefore, we have confidence 
that deviations from spherical symmetry in HD 153919 
and/or the effects of X-ray irradiation are negligible for 
the purposes of spectroscopic modeling. 

3. Mass determination for the system components 

Since no X-ray pulsations have been convincingly mea- 
sured for 4U1700-37 - and hence no determination of 
a^^sini is possible - the orbital solution cannot be uniquely 



determined. However, following Heap & Corcoran (1992) 
and Rubin et al. ( 1996] ) we may estimate the mass of the 
companion using a Monte Carlo method. The mass of the 
companion can be calculated using a series of equations 
relating eclipse and orbital parameters. 

The Roche lobe filling factor il is defined by 
= VLRl (1) 
where i?* is the radius of the 0-star and Rl is the Roche 
lobe radius which is related to the semi-major axis of the 
system a and the mass ratio q (= Mast/Ma,) by 



Rl 



= A + B\ogq + C{\Qgqf 



(2) 



and the coefficients A, B and C are 



0.398 - 0.026r^ 



0.004r3 

o.oi5r3 



A 

B = -0.264 + 0.052r2 
C = -0.023 - 0.005r2 

where T is the ratio of the rotational angular frequency 
of the companion to its orbital angular frequency (Rap- 
paport & Joss 1983 ). 

The radius of the 0-star is related to the semimajor 
axis by the inclination i and eclipse semiangle by 

— = \J cos^i + sin^icos^^E (3) 
while the companion mass function / is given by 



and / is related to orbital parameters by 
/ = 1.038 X IQ-'^K^P{1 - 6^)3/2 



(4) 
(5) 



where is the radial velocity semi-amplitude in km s^^, 
P the period in days and e the ellipticity. Finally by Ke- 
pler's third law 



We combine these equations in a similar way to Rubin 
et al. ( |1996| ) to obtain 



Rl^ 



Rla 



17.81P^/^f/^{l + qfcos^0E 



(7) 



75.19(1 -t-g)M^F^ 



(6) 



where Rlu = Rl/ci. This equation can then be solved 
numerically for q. 

Values of various system parameters (listed in Table |2|) 
are selected randomly either from a gaussian distribution 
(if observed) or uniformly if constrained between certain 
values. Equation |^ is solved for q which gives a from 
Eq. H which then gives Ma, from Eq. ^. Consistency can 
be checked by requiring sin z < 1 and i > 55 degrees (e.g. 
Rubin et al. |1996D . 

This procedure is followed 10^ times to gain a distribu- 
tion of Mx and for a considerable number of possible 
parameter combinations (noting that as expected there is 
a very strong positive correlation between the two masses) . 
We find that M^; = 2.44 ± 0.27 Mq. As shown by the his- 
togram of Ma; in Fig. 4 the distribution is very asymmet- 
ric beyond the la limits (determined by the 16*'' and 84*^* 
percentile of the cumulative distribution function), with 
only 3.5 per cent of the sample having a mass of less than 
2 Mq, and none less than 1.65M0 (which is significantly 
higher than the upper limit to the ra nge f ound for binary 
pulsars by Thorsett & Chakrabarty 1999). We note that 
none of the 10^ trials were rejected from inclination con- 
straints suggesting that the range of stellar radii adopted 
for the modeling are unlikely to be significantly in error 
(which, for the fixed eclipse length, would lead to unphys- 
ical solutions for the orbital inclination). 

The errors on Ma, are significantly smaller than pre- 
vious work (e.g. Rubin et al. 1996 ) due to the far more 
stringent limits on R^, which constrain the orbital and 
eclipse parameters far more strongly. This is not surpris- 
ing as the eclipse parameters are used to work out the or- 
bital parameters and the eclipse constraints rely strongly 
on the 0-star radius. 

Fig. 5 shows the 0-star mass distributions around 
M* = 58±11 Mq. Again the distribution is anti-symmetric 
with 32 per cent of trials between 50-60 Mq, 26 per cent 
between 40-50 Mq and only 2 per cent less than 40 Mq. 
Therefore, the mass implied for HD 153919 appears to be 
consistent with both that expected from its spectral classi- 
fication and relevant evolutionary tracks (see Fig. 6), and 
that suggested by its high terminal wind velocity (Sect. 
2). Additionally the \ogg determined from the Hei and 
Balmer line wings (Sect. 2.2) indicates a minimum mass 
of 50 Mq (and maximum of ^60 Mq), again fully con- 
sistent with the results of the Monte Carlo simulation. 
Therefore, given the consistency between mass estimates 
based on spectral type, evolutionary tracks (when com- 
pared to the stellar temperature and luminosity derived 
from modeling) , surface gravity and the Monte Carlo sim- 
ulations, we have confidence that the mass of HD 153919 
lies in the range 50-60 Mq . This resolves the problem that 
the star is undermassive by a factor of ~2. 
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However, the mass of the compact companion is more 
problematic given that it is significantly in excess of the 
observed mass range for NS, but apparently considerably 
lower than those found for BH candidates (e.g. Fig. 7). If 
a minimum mass of 50 M0 is adopted for HD 153919 the 
minimum value of that may be obtained is 1.83 M0, 
while for values of Mq between 50-60 Mq only 0.17 per 
cent of trials result in M^; <2 M©. This will be returned 
to in Sect. 5. 

Recent reanalysis of spectroscopic data by 
Hammerschlag-Hensberge et al. (in prep.) suggests 
that the eccentricity of the orbit is somewhat uncertain, 
and that the orbital velocity curve is equally well fit 
by an orbit of eccentricity e ~ 0.22 ± 0.04 as it is by a 
circular orbit. In order to address this uncertainty we 
modified the above equations for the more general case 
of an elliptical orbit and repeated the simulations with 
e — 0.22 ± 0.04. This resulted in significantly higher 
masses for both components, with M*=70±7 M© and 
M^=2.53±0.2 Mq. Therefore, the mass of the O star in 
the case of an elliptical orbit is significantly higher than 
expected for an 06.5 Iaf+ star (only 0.002 per cent of 
the trials result in a mass <50 M©, and 5 per cent give 
a mass between 50-60 M©). Such high values for are 
inconsistent with the measured logg and we note that 
95 per cent of trials are rejected due to the inclination 
constraints, suggesting that a low eccentricity solution is 
more likely. 

If such extreme values for are adopted, the mass of 
the compact object is still less than that observed for the 
lowest mass black hole candidate known (~4.4 M©; Sect. 
5) and remains significantly greater than any known neu- 
tron star. Indeed, the lowest mass estimates for both com- 
ponents were derived in the case of a circular orbit; there- 
fore the value of Mj;=2.44±0.27 M© represents a lower 
limit for the mass of the compact object^, and we suggest 
that these results favour a low eccentricity solution for the 
orbit (we note that the orbital eccentricity of Vela X-1 is 
overestimated from optical observations when compared 
to th e value derived from timing analysis, cf. Barziv et al. 
2OO1I) . 



4. Evolutionary history of 4U1700-37 

The stellar parameters for both primary and compact 
companion determined via non-LTE modeling and Monte- 
Carlo simulation raise many important questions regard- 
ing the evolution of single and binary massive O stars and 
their ultimate post-SN fate. However, such questions are 
complicated by the uncertain evolution of massive stars 
after leaving the main sequence and the role that binarity 
and associated - possibly non-conservative - mass transfer 



^ Values of Ma; <2 Mq are only obtained when the radial 
velocity semi-amplitude, Ko, is >2 a below the observed value, 
to obtain M^: < 1.6 M0 would require Ko to be wrong by several 
sigma. 



2.5 







Fig. 4. Histogram of the results of the Monte Carlo simu- 
lations for the mass of the compact object in 4U1700-37. 
The results indicate a mass in the range of 2.44±0.27 i?© 
with only 3.5 per cent of simulations indicating masses of 
less than 2 M©, and none <1.65 M©. 




100 



Fig. 5. Histogram of the results of the Monte Carlo simu- 
lations for the mass of the 06.5Iaf-f primary HD 153919 
- the results indicate a mass in the range of 58±11 i?© 
consistent with evolutionary predictions and the mass es- 
timated from our determination of log (7=3.45-3.55 (Sect. 
2.2). 

plays in modifying this; for instance the time at which the 
hydrogen rich outer layers are lost exposing the helium 
core plays a critical role in determin ing t he final pre-SN 



mass of the star (e.g. Brown et al. ^2001 and references 

therein). 

Heap & Corcoran ( 1992 ) propose an initial 
80 Mq-I- 40 M© binary system with subsequent evolution 
via case B mass transfer^ after 2.6x10^ yrs, proceeding 
for 10^ yrs (via Roche-lobe overflow; RLOF). After the 

We adopt the nomenclature used in Wellstein & Langer 



(1999 and references therein), with Case A, B & C evolution 



J. S. Clark et al.: Physical parameters of 4U 1700-37 



9 




Fig. 6 . The location of HD 153919 in the Hertzspmng- Russell diagram. The evolutionary tracks of Lejeune & Schaerer 
( 200C ) for stars of a given initial mass are indicated (up to the phase of core-helium burning). The numbers along 
the tracks show the decrease in M* with time due to wind losses. The plus signs indicate the end of core-hydrogen 
burning. The diagonal dashed lines are lines of constant radius. 



Parameter 



Value 



F 
ft 

e 

9b 
P 



21.4 - 23.2 Rq 
0.5 - 1.0" 
0.8 - 1.0" 
0.0'' 

28.6 ±2.1 degrees" 
3.411581 ± 2.7 X 10"^ days" 
20.6 ± 1.0 km S-" 



Table 2. Physical parameters of 4U 1700-37. Those with 
a ± have a gaussian error distribution while those without 
are assumed to have a uniform distribution. ° Rubin et al 
( |1996D , ''Hammerschlag-Hensberge et al. (in prep.) 



mass transfer the initially more massive star has lost 
enough material due to the combination of wind driven 
mass loss and the brief period of non-conservative RLOF 



corresponding to mass transfer during core hydrogen burning, 
after core hydrogen burning but before core helium exhaustion, 
and after core helium burning, respectively. 



to become a WR star, which subsequently explodes as a 
SN. 



Based on their identification of the Sco OBI asso- 
ciation as the birthplace of HD 153919, Ankay et al. 
(2001) propose a lower initial mass of the SN progenitor 



(>30_2o ^^o) based on the turnoff mass for the proposed 
6±2 Myr age of Sco OBI at the time of the supernova. As- 
suming conservative Case B mass transfer they derive an 
initial mass of at least 25 M0 and suggest the short orbital 
period Wolf-Rayet (WR) binary CQ Cep/HD 214419 as a 
possible example of the progenitor system. However they 
note that the assumption of conservative mass loss might 
be incorrect and highlight the non-conservative scenario of 



Wellstein & Langer (1999). Such a scenario is attractive 



since the loss of substantial quantities of mass and an- 
gular momentum naturally lead to short period binaries 
(assuming both components do not merge). However, com- 
mon envelope evolution is poorly understood and therefore 
somewhat limits our ability to quantitatively reconstruct 
the pre-SN evolution of the binary. 
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Despite these uncertainties we can address the general 
evolution of the binary in some detail. Our present mass 
estimate for HD 153919 of = 58^11 M© suggests a 
mass for the SN progenitor of the order of <;60 M0, at the 
upper range of that proposed by Ankay et al. ( ^001 )^. The 
short orbital period of HD 153919/4U1700-37 favours non 
conservative evolution probably via case B mass transfer. 
While mass loss via the stellar wind of the SN-precursor 
during the WR phase will lead to a widening of the orbit, 
a favourable SN kick may overcome these problems. 

The alternative Case C evolution appears unlikely 
given that the SN will occur several thousand years after 
the end of mass transfer/loss. This time period appears 
unrealistically short given that this will not allow suffi- 
cient time for the helium mantle to be removed to expose 
the C/0 core (i.e. the star will not pass though a WC 
phase). Therefore the SN progenitor would have a large 
mass at the point of SN; for a SN progenitor with an in- 
tial mass of 60 M© we might expect the mass to be of 
the order of 30 Mq (the maxium helium core mass) . This 
implies the loss of a very large quantity of material in the 
SN, contradicting the estimate of Ankay et al. (2001) that 
only ~9 M© was lost in the SN (and we might also ex- 
pect such a scenario to lead to a very high mass compact 



object, cf. Brown et al. 2001). 

The high mass implied for the SN precursor suggests 
that such a star would be likely to evolve through a Lu- 
minous Blue Variable (LBV) phase rather than a Red Su- 
pergiant (RSG) phase on its way to becoming a WR star 
(stars with initial masses ^40 M© likely avoid the RSG 
phase) . Such an evolutionary path is likely to prevent mass 
trans fer onto HD 153919 via RLOF (Wellstein & Langer 
1999| ). Significant accretion of material by HD 153919 via 



RLOF seems implausible in any case, as this would lead to 
a larg e orbital separation and period (Wellstein & Langer 
199£). Instead, the formation and subsequent ejection of a 



common envelope (despite the SN precursor avoiding the 
RSG phase) and simultaneous reduction in orbital period 
and binary separation is suggested. Such a scenario there- 
fore implies that the present day mass of HD 153919 forms 
a lower limit to the mass of the SN precursor, subject to 
the possible accretion of a small quantity of material di- 
rectly from the wind of the SN precursor (see below). 

After the ejection of the outer layers of the SN pre- 
cursor we are left with a short orbital period WR+0 star 
binary. Support for this scenario is provided by the possi- 
ble overabundance of carbon (or rather the lack of signifi- 
cant C depletion as might be expected for CNO processed 
material; Sect 2.2) in HD 153919. The carbon rich ma- 
terial must have originated in the SN precursor during a 



^ We note however, that the difference in the life time of 
a 60 M0 star (3.8 Myr) and - for example - a 120 Mq star 
(3.0 Myr) is smaller (0.8 Myr) than the uncertainty in the life 
time of a 60 Mq star (1 Myr), suggesting that higher progenitor 
masses may be possible. 



carbon rich WC stage, independently suggesting a rather 
high initial mass for the SN precursor. Subsequent mass 
transfer would then have to occur via direct wind fed ac- 
cretion, with the wind of the WC star impacting directly 
on the surface of the O star. Despite the high mass loss 
rate of the O star (M=9.5xl0-6 M©yr-i; Table 2) the 
possibility of such accretion is suggested by hydrodynam- 
ical simulations of colliding wind binaries (Gayley, Owocki 



Cranmer 1997; Dessart, Petrovic & Langer, in prep.). 



We may exclude the overabundance in nitrogen orig- 
inating via direct wind accretion during the WN phase 
of the SN precursor. For a nitrogen overabundance in 
HD 153919 of ^9, the nitrogen mass fraction, Xjv=0.01, is 
the same value as is found in the winds of WN stars which 
implies that HD 153919 would have to accrete 18 M© from 
the SN precursor during the WN phase to produce such an 
overabundance. Given that this is unreasonably high, we 
suggest that the excess nitrogen probably originated from 
internal, rotational mixing (nitrogen overabundances are 
not unusual for O stars). 

Several short period WR-l-0 star binaries are known 
and could provide analogues to the precursor of the 
present binary configuration. At present CQ Cep does 
not fit particularly well (MvKAr=21 M©, Mo9=26 M©, 
P=1.6 days) - in a few 10^ yrs when the WN star 
has lost mass and has evolved into a lower mass WC 
star (and the period has lengthened) it may provide 
a better model for the system. However, of the known 
WC+0 binaries, HD 63099 {Mwc^9 M©, Mo7=32 M©, 
P=14 days) could evolve into a HD 153919/4U1700-37 
like binary if the SN kick is favourable. Other known 
WR-l-0 star binaries that could form similar systems are 
HD 152270 (MvKC=llM©, Mo5-8=29 M©, P=8.9 days) 
and HD 97152 (Mvkc = 14 M©, Mo7=23 M©, P=7.9 days); 
in all three cases the present mass of the WC star is >9 M© 
as suggested for the SN precursor by Ankay et al. (2001) 
on the basis of the present mass of the compact companion 
and the current space velocity of the system. Therefore, of 
the si x WC-l-0 binaries with known masses (van der Hucht 
^OOlD thr ee are found to have parameters consistent with 
the presumed pre-SN stage of 4U 1700-37. 

The results of the Monte-Carlo simulations suggest 
that the SN formed a compact object (see Sect. 5) with a 
mass in the range 2.44±0.27 M©. As massive binary mod- 
els show that Case B primary stars of larger initial mass 
evolve to a larger final mass, and that Case A primaries 
end u p less massive than Case B's (cf. W ellstein & Langer 
1999| ; Wellstein, Langer & Braun ^00l| ) we can conclude 



that the remnants of all Case A/B primaries initially less 
massive than about 60 M© are less massive than about 
2.5 M©. 

Theoretical predictions suggest that a 60 M© star in 
a close binary system is capable of producing either a low 
or a high mass compact object depending sensitively on 
the wind mass loss rate adopted for such a star during 
its WR phase; a variation of only a factor of three in the 



J. S. Clark et al.: Physical parameters of 4U 1700-37 



11 



Neutron Stars 

1.4M.. 



Vela X-1 

Her X- 1 

4U1538-5a 

Cen X-3 

LMC X-4 

SMC X-1 

J0045-73 

B1855+09 

B180a-07 

J1713+07 

,11012-1 53 

B159B+49(C) 

B1598+49(P) 

BZ303+46{C) 

B2303+46(P) 

Bai27+U{C) 

B2127+11{P) 

B1534+12{C) 

B1534+12{P) 

B1913+16{C) 

B1913+16{P) 



Black Holes 



lOM.. 



H4U1700 37 



LMC X-1 
LMC X-3 
^Cyg X-1 
J1550-564 
4U1543-47 
V4641 Sgr 
J165S-40 
J0422+32 
A0620 00 
K Hus 91 
N Oph 77 
02000+25 
V404 Cyg 

jiiia-4ao 

CHS1915+ 105 
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Fig. 7. Mass d istribu tion for neutron stars and black holes 
(after Charle s 1998 ). Neutro n star masses are from Ash 
et al. (1999), Barziv et al. ( 2001 ), va n Ke rkwijk et al. 



( |1995D and Thorsett & Chakrabarty ( [l999| ). Black hole 
masses provided by Charles (priv. comm.). Error bars for 
systems other than 4U1700-37 are 1 sigma errors; see Sect. 
3 for a discussion of the errors associated with the mass 
of 4U1700-37 but note that the probabihty of M<2Mo is 
less than 3.5 per cent, and no trials results in M<1.65M0. 



WR mass loss rate leading to compact object masses in 
between 1.2 and 10 Mq (Fryer et al. 



2002). This result 



shows that given the present uncertainties in WR mass 
loss rates, the relatively low mass found for the compact 
star in 4U1700-37 is not in conflict with the evolutionary 
scenario proposed above (and also argues for a relatively 
high WR mass-loss rate). 

Therefore, given the low mass of the compact com- 
panion in 4U1700-37, it seems to be difficult to explain 
any of the (high mass) galactic black hole binaries as be- 
ing produced through the Ca se A/ B channel (cf., Porte- 
gies Zwart, Verbunt & Ergma 1997) except for those with 
the most massive SN-progenitors (Brown et al. |2001 ). In- 
deed, evolutionary scenarios invoking Case C mass trans- 
fer (Brown, Lee, & Bethe |l999| ) seem to be required to 
explain the high-mass black holes in low-mass X-ray bina- 
ries. 



5. The compact companion 

At present mass determinations exist for 36 compact ob- 
jects, of which 21 are neutron stars and the remainder 
black hole candidates (Fig. 7). Thorsett & Chakrabarthy 
( |1999D found that the masses of neutron stars are clus- 



tered in a remarkably narrow range (mean of 1.35Mq and 
a standard deviation of 0.04 Mq). However, recent anal- 
ysis of Vela X-1 by Barziv et al. (|200l[ ) suggests that the 



pulsar has a mass of l.STtolr M0I], while Orosz & Kuulk- 
ers ( |1999D find a mass of 1.78±0.23 Mq fo r Cygnus X-2 



(however Titarchuk & Shaposhnikov 2002 have recently 
proposed M=1.44±0.06 M© on the basis of Type-I X-ray 
bursts). 

Based on assumptions about the origin of kilohertz 
quasi-periodic oscillations Zhang et al. (1997) suggest that 
several LMXBs may also contain massive (~2Mq) neutron 
stars, resulting from the accretion of substantial amounts 
of material over long (10^ yrs) periods of time. However 
their putative descendants, radio pulsar -|-white dwarf bi- 
naries, provide no evidence for massive neutron stars ( 
Barziv et al. 2001). The same is true for the Be/X-ray bi- 
naries, which have evolved from lower mass systems than 
the OB-supergiant HMXBs. Only in the latter systems ev- 
idence has been found for massive neutron stars and black 
hole candidates. 

Masses for a number of black hole candidates have also 
been determined; lower limits to their masses comfortably 
exceed 3M0. Indeed several objects appear to have masses 
>10Mo (e.g. 14±4Mo for GRS 1915+105; Greiner, Cuby 
& McCaughrean |200l| ), with most typically ^6Mq. At 
present the binary system with the lowest mass candidate 
black hole is Nova Vel (4.4 Mq; 



199c) 



Fillipenko et al. 
Given that the mass of the compact object in 4U1700- 
37 lies outside the present observational range for both 
neutron stars and candidate black holes, is it possible to 
increase/decrease the mass of the compact object such 
that it is consistent with cither type of object? For the 
case of a circular orbit we find no solutions consistent with 
Ma; ^ 4.4 Mq. Very eccentric orbital solutions do allow 
values of M^: in this range although we note there is no 
physical motivation for them. Such solutions imply values 
of Mo greatly in excess of 100 M0, well above current 
evolutionary predictions for the mass of an 06.5 Iaf-|- star. 
Stars with masses in excess of 100 Mq are instead expected 
to evolve to H depleted WR stars via a H rich pseudo-WNL 
phase where their powerful stellar winds mimic the spectra 
of more (chemically) evolved stars of lower masses. 

Such a large mass would also be inconsistent with the 
constraints imposed by the surface gravity (Sect. 2.2) and 
the relationship between stellar mass and the terminal 
wind velocity. Finally the result would imply that even 
very massive stars (given that the initial mass of the SN 
progenitor had to be significantly in excess of the present 
mass of HD 153919) leave relatively low mass remnants 
post supernova, presenting significant problems for the ori- 
gin of heavy (>10 Mq) black holes such as e.g. Cyg X-1. 



® Systematic excursions in the radial velocity curve for this 
system complicate this determination and prevent an unam- 
biguous confir matio n of the mass estimate for the neutron star 
(Barziv et al. |200l|) 
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Therefore we conclude that M^; appears to be inconsistent 
with the range of masses of known black hole candidates. 

Given the stringent lower limits derived in Sect. 3 it 
also appears difficult to bring the mass of the compact ob- 
ject into line with the range of masses found by Thorsett 
& Chakrabarthy ( 199£ ). Inspection of the evolutionary 
tracks in Fig. 6 suggest that stars with initial masses of 
the order of 60 Mq initially evolve redwards before re- 
turning bluewards after undergoing significant mass loss, 
most likely during an LBV phase. While we note that the 
behaviour of stars in this short lived phase is very uncer- 
tain it is unlikely that we could be observing HD 153919 
after such an excursion, since we would expect significant 
chemical enrichment - the H rich mantle having been lost 
- which is not observed. Likewise, the mass constraint im- 
posed by the determination of the surface gravity also ap- 
pears to exclude this scenario. 

Furthermore, such a low value for Mq and would 
reintroduce the problem of HD 153919 being undermassive 
for its spectral type by a factor ^2. While the primaries in 
some HMXB systems are found to be undermassive (e.g. 
Gen X-3 and LMG X-4; Kaper ^OOlD this is attributed to 
mass loss via RLOF - wind fed systems do not show this 



effect (Kaper |200lD . Given that 4U1700-37 is currently a 
wind fed system and is yet to evolve into a RLOF system 
this could not explain such a mass discrepancy. Indeed, 
given the evolutionary constraints imposed by the present 
orbital period it is likely that no significant mass transfer 
has occured onto HD 153919 during its lifetime, and it will 
have evolved as if it were an isolated star. 



5.1. Implications of an intermediate mass 

If we accept M2;=2.44±0.27 Mq - as implied by the sim- 
ulations - the nature of the compact companion remains 
uncertain. Conflicting claims as to the nature of the ob- 
ject have been made on the basis of the X-ray spectrum 
of 4U1700-37 (Sect. 1). While we cannot discriminate be- 
tween the twin possibilities of massive neutron star or low 
mass black hole from our present measurements we note 
that consideration of the masses of both components of the 
binary appear to exclude the possibility that stars with 
masses of '~60 Mq can produce 5-10 M0 black holes via 
case A or B evolution. 

If the compact object in 4U1700-37 is a black hole it 
confirms Brown et al.'s prediction of the existence of 'low 
mass black holes' (based on their "soft" equation of state, 
Brown et al. 1996|) , while if the object is a neutron star the 
high mass would severely constrain the equation of state 
of matter at supra-nuclear densities. 

The relationship between the mass of a neutron star 
and its central density is calculated by integration of 
the Tolman-Oppenhe imer-V olkoff (TOV) equation (Op- 
penheimer & Volkoff 1939) which is the relativistic ex- 
pression for hydrostatic equilibrium. In order to perform 



the integration it is neccesary to understand the equation 
of state of the degenerate nuclear matter in the star. 

Because of their non-perturbative nature, strong inter- 
actions between nuclei are extremely difficult to calculate 
even under normal conditions. However, there are several 
models of the internuclear potential in the literature which 
have achieved much success in modelling nuclei e.g. Stoks 
et al. ( p^94|) , Whinga et al. ( |l995| ) and Machleidt et al. 
( |l996| ) . One of the most successful and up to date of these 
models has been applie d to t he neutron star equation of 
state by Akmal et al. ( |1998| ) yielding a maximum mass 
of between 2.2 and 2.4M0 for a neutron star made com- 
pletely of normal nuclear matter. 

There is also the possibility of a QGD phase transition 
occuring in the centre of the neutron star. The large chem- 
ical potential has a similar effect to a large temperature on 
the QGD coupling constant. Gonsequently the interquark 
coupling can be reduced to the point where deconfinement 
occurs and nuclei dissolve into quark matter. The presence 
of quark matter has the effect of softening the equation of 
state which leads to a lower possible maximum mass for 
the neutron star. The energy scale at which deconfinement 
occurs can be parameterised by the QGD bag constant B, 
a phenomenological parameter representing the difference 
in energy density between the vacua of hadronic and quark 
matter {B « 120 - 200 MeV fm'^). Inclusion of a phase 
transition to such a mixed state reduces the maximum 
mass of the neutron star to ~2 Mq for B — 200 MeV 
fm''^ and - 1.9 Mq for B = 122 MeV fm'^ (Akmal et al. 
19981 )11. 

If the neutron is rotating rapidly, the TOV equation 
ceases to be a valid approximation and one must drop the 
assumption of spherical symmetry in the metric. The ef- 
fect of the rotation will increase the allowed mass of the 
neutron star as one might expect, but it is shown (Heisel- 
berg & Hjorth- Jensen 1999 ) that even with rotation one 
cannot obtain a neutron star with a mass much higher 
than those listed above without using an equation of state 
so stiff that the sound speed in the neutron-star interior 
becomes superluminal. 

To summarize this subsection, state of the art models 
for nuclear equations of state which include the effects of 
three nucleon interactions marginally allow the existence 
of a 2.4 Mq neutron star. However, the existence of such a 
star would place severe constraints upon the onset of new 
physics at high hadronic densities. For instance, in the 
model where a quark matter core is expected to develop, 
the bag constant B would have to be considerably larger 
than 200 MeV fm^"^ in order for such a star to be viable. 



^ Hyperons such as A, E and A may also be produced at 
high density and it is thought that their presence will also 
soften the equation of state and reduce the ove rall p ermitted 
neutron star mass (Heiselberg & Hjorth- Jensen 1999). 
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6. Summary 

We have performed a sophisticated NLTE analysis on 
the 06.5 Iaf+ star HD 153919, the primary of the 
HMXB 4U1700-37 and used the resuhs to constrain 
the masses of both components of the system via 
a Monte Carlo simulation. Our NLTE model atmo- 
sphere analysis leads to parameters for HD 153919 
of log(LjL0)=5.82±O.O7, Teff=35,000±1000 K, 
R,=21.9lo;^ R-0> M=9.5xl0~s Mq yr"!, and an 
overabundance of nitrogen and possibly carbon over solar 
metallicities. Combined with the short orbital period 
of the system this implies a common envelope phase of 
pre-SN evolution - despite the mass of the SN progenitor 
apparently precluding a RSG phase - leading to the 
formation of a close WC-l-0 star binary, with the carbon 
enrichment of HD 153919 a result of the impact of the 
stellar wind of the WC star on the surface of the O star. 

The Monte Carlo simulations result in masses for 
the O star and compact object of M*=58±ll M© and 
M;^,=2.44±0.27 Mq, with none of the 10^ trials result- 
ing in Ma; <1.65 Mq, while only 3.5 per cent of the tri- 
als result in M^, <2 Mq. Given that no significant mass 
transfer via RLOF has occured this implies that the ini- 
tial mass of the SN precursor must have been ^60 Mq. 
Thus even very massive stars can effectively 'melt down' 
to leave rather low mass post-SN remnants. Equally, the 
masses of both components imply that it is impossible for 
stars of ~60 Mq to leave 5-10 Mq remnants via Case A or 
B evolution, suggesting that most high mass black holes 
are instead formed via Case C mass transfer. 

The mass of the compact object is found to lie in be- 
tween the range of masses observed for neutron stars and 
black holes. Given that M* and M^; are strongly corre- 
lated, forcing consistency between Ma; and either type of 
object results in significant discrepancies between M^,, and 
evolutionary predictions for the mass of HD 153919. 

In order to produce consistency with the range of 
masses observed for neutron stars the O star has to be sig- 
nificantly undermassive for its spectral type and luminos- 
ity class. While theoretical evolutionary tracks for massive 
(^60 Mq) stars suggest that after a redwards excursion 
the star will evolve bluewards again with a substantially 
reduced mass, such a star would show significant chemi- 
cal enrichment (and H depletion) which is not observed. 
Equally, the surface gravity determined from modeling ex- 
cludes such an anomolously low mass. 

Forcing consistency between M^, and the masses of 
known black hole candidates {Mx ^ 4.4 Mq) results in 
M* ^100 Mq. Stars of such extreme masses are not ex- 
pected to go through an O supergiant phase, rather evolv- 
ing into H depleted WR stars via a H-rich pseudo WR 
phase, where their high mass loss rate simulates the spec- 
trum of a WR star. Once again, the constraint implied by 
the surface gravity also appears to exclude this possibility. 



We are therefore left with the conclusion that no solu- 
tion is fully consistent with present expectations for stellar 
evolution and the chemical abundances and surface grav- 
ity of HD 153919. If the compact object has a mass con- 
sistent with the observed range of neutron star masses, 
the O star is significantly undermassive, while if it consis- 
tent with the lower limit to black hole masses the O star 
is overmassive by a similar (or larger) factor. Finally if - 
as the Monte Carlo analysis implies - the probable mass 
of the O star is consistent with evolutionary predictions 
and the measured surface gravity, the mass of the compact 
object lies in between the two alternatives. 

While our results do not allow us to distinguish be- 
tween a massive neutron star or a low mass black hole, 
the existence of a neutron star of mass in the mass range 
2.44±0.27Mq would significantly constrain the high den- 
sity nuclear equation of state and provide details about 
the QCD phase transition complementary to information 
about the temperature induced transition which will be 
obtained at RHIC and LHC. Phenomena which might oc- 
cur deep in the star such as the appearance of hyperons 
or a quark matter core would be strongly constrained as 
the existence of these phases might result in an equation 
of state too soft to support such a high mass star. 
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